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ANALYSIS OF STRESSES AT THE BORE OF A DRILLED BALL 
OPERATING IN A HIGH-SPEED BEARING 
by Harold H, Coe and John E. Lynch 
Lewis Research Center 

SUMMARY 

A computational analysis was performed on two drilled- ball configurations to deter- 
mine the effect of a change in geometry on the three-dimensional stress distributions in 
the balls. A finite- element technique was used to calculate the stresses in a regular 
drilled ball (50 percent weight reduction) and in a drilled ball with a stiffening web 
(46 percent weight reduction). The drilled ball was 20. 6 millimeters (0. 8125 in. ) in di- 
ameter and had a 12. 6 millimeter- (0. 496-in. -) diameter concentric hole. The stiffening 
web used was 1.5 millimeters (0,06 in.) thick. 

The calculations were made for conditions of a drilled ball operating in a 125- milli- 
meter- diameter -bore bearing at a shaft speed of 24 000 rpm (a DN value of 3 million) and 
at thrust loads of 11 000 and 22 000 newtons (2500 and 5000 lb). These were the approx- 
imate conditions under which an experimental drilled-ball bearing of the same size ex- 
perienced a ball failure due to flexural fatigue. 

The analysis shows that a large reversing stress occurs in the bore of the ball as the 

ball rotates through one complete revolution. At the 22 000 newton- (5 000- lb) load con- 

£* 

dition, the maximum tangential stress at the bore ranges from 470x10 newtons per square 
meter (70 000 psi) in tension to about 150X10^ newtons per square meter (20 000 psi) in 
compression. 

The addition of the stiffening web reduces the change in tangential stress to about one- 
half that in the regular drilled ball or 300x10'^ newtons per square meter (45 000 psi) at 
the 22 000- newton- (5 000- lb-) thrust load condition. 


INTRODUCTION 

Recent trends in gas turbine design have resulted in a requirement for higher shaft 
speeds and larger shaft diameters (ref. 1). While bearings in current aircraft turbine 



engines operate in the range 1.5 to 2.0 million DN (bearing bore in mm multiplied by 
shaft speed in rpm), bearings in future engines may operate at DN values of 3. 0 million 
or higher. When ball bearings are operated at high DN values, however, centrifugal 
forces on the outer race produced by the balls can become significant. The resulting 
increase in Hertz stress at the outer-race contacts can seriously shorten bearing 
fatigue life (ref, 2). 

A possible solution to this high-speed problem is to reduce the mass of the ball and 
thereby reduce the centrifugal forces (ref. 3). Therefore, bearii^s with spherically 
hollow or cylindrically hollow (drilled) balls have been tested (refs. 4 to 8). In these 
tests, the 17. 5-millimeter- (0. 6875-in. -) diameter spherically hollow balls with a 1.5- 
millimeter (0, 060-in.) wall thickness (refs, 4 and 5) failed by flexural fatigue in the 
area of the electron-beam weld. The drilled-ball bearings performed satisfactorily at 
DN values up to 3 million, as described in references 6 and 7. After accumulating some 
running time, however, the bearings experienced ball failures, as reported in refer- 
ence 8. It was concluded that these ball failures were caused by cracks that, because 
of flexural fatigue, initiated in the bore. To understand these drilled-ball fatigue fail- 
ures better and to evaluate a proposed ball design change, three-dimensional stress 
distributions were required. 

Therefore, this report describes the computation and analysis of the stresses for a 
drilled ball and for a drilled ball with a stiffening web, A finite- element technique was 
utilized for the computation. While the stresses were computed for the complete ball, 
the analysis is mainly concerned with the stresses at the bore. 

The analysis was made for a 20. 6-millimeter- (0. 8125-in. -) diameter drilled ball 
under the conditions of operation in a 125-millimeter-diameter-bore ball bearing with 
an inner-race speed of 24 000 rpm and bearing thrust loads of 11 000 and 22 000 newtons 
(2500 and 5000 lb). The drilled ball had a 12. 6- millimeter- (0. 496-in. -) diameter hole 
(50-percent weight reduction). The drilled ball with a stiffening web had two 12. 6-milli- 
meter- (0. 496-in. -) diameter holes and a central web, 1. 5- millimeters (0, 06 in. ) thick, 
between the holes in the center of the ball. The weight reduction for this ball was 46 per- 
cent. 


ANALYTICAL PROCEDURE 

It was decided to calculate the stresses that exist in a ball operating under the con- 
ditions noted in reference 8 prior to failure. The 125- millimeter- bore bearing used 
in reference 8 and considered in this report is shown in figure 1. The specifications for 
this bearing are shown in table I, The bearing operating conditions as well as the cal- 
culated values for ball loads, contact armies, and speeds are shown in table II, 
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Also considered for analysis was a drilled ball with a stiffening web. The condi- 
tions calculated for a bearing with the webbed ball are also shown in table II, Cross 
sections of the drilled ball and of the ball with a stiffening web are shown in figure 2. 

The web thickness used was 1. 5 millimeters (0. 06 in. ). 

Because of the centrifugal force generated by the balls orbiting about the bearing 
axis, the inner- and outer-race contact loads and angles become dissimilar, as shown 
in figure 3. It was for this load configuration that the stresses in the balls were deter- 
mined. 

The spherical balls contain a cylindrical hole. A cylindrical coordinate system is 
used to locate points within the ball, with the z-axis collinear with the hole centerline 
(fig. 3). 

In this coordinate system, a positive radius vector is directed away from the z-axis 
and is normal to it. Although the ball geometry is axisymmetric about the z-axis, the 
loads, deflections, and stresses are dependent on the angle 0, which is measured as 
shown in figure 3. 


STRESS ANALYSIS 

A finite- element technique was utilized to calculate the stresses in the ball. The 
ball was divided into axisymmetric finite rir^ elements as shown in cross section in 
figures 4 and 5. 

The stress calculations were performed on a digital computer usii^ a ring finite- 
element program. This program analyzes axisymmetric solids of revolution with ar- 
bitrary shape which are subjected to unsymmetrical loadii^s. It determines elastic 
deformations, stresses, and strains in a three-dimensional solid which is modeled as 
an assembles of rings of regular cross section. 

All loads are approximated by Fourier series load terms in the tangential coordi- 
nate d. The three-dimensional equations of elasticity are uncoupled into sets of two- 
dimensional equations, one set of equations for each Fourier series loading. 

Each set of equations is solved for the Fourier nodal deflections and element stress 
components for that Fourier loading. These Fourier stress components are finally 
combined to give the r, z, and $ directional stresses for each element at any selected 
tai^ential angle. 

The number of Fourier terms is arbitrary, but is restricted to 25 for this investi- 
gation. A greater number of terms would improve the angular representation of con- 
tact loads and contact stresses, at the expense of increased computer time, with little 
improvement in representing bore stresses according to Saint Venant's principle. 
Twenty-five Fourier terms appears to be a practical compromise. 
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Contact Force Representation 


The load at a ball contact ’’point' ♦ is actually distributed over an area which depends 
upon the load magnitude. If treated as a pressure, the load is a function of r, z, and Q 
within the contact area. However, in the present analysis, the load is represented as a 
discrete function of r and z and as a continuous function of B in the following manner. 
The load is first distributed over a small number of nodal points located in the r, z-plane; 
these nodes are adjacent to the contact "point." The r or z load component at each 
nodal point (0) is then approximated as a delta function of angle 



(0) = F.6(0 - p) “ lim F. 


sin c(0 - p) 
v(0 - p) 


( 1 ) 


where i = r for the radial force component at the node, i = z for the axial force component 
at the node, p = 0° for outer-race loads, and p = 180° for inner-race loads. This ai^- 
ular variation is shown in figure 6. 

A finite value of c was selected to define the angular "width" of the contact area, 
which was then reasonably close to a normal contact dimension. Since the numerical 
method used for the stress analysis utilizes a Fourier coefficient approach, the Fourier 
load coefficients must also be obtained. The load component at each nodal point is rep- 
resented in a Fourier cosine series as 




n=24 

/ /n ^ T 

n=l 


( 2 ) 


where the n*^^ Fourier load coefficient is 


f+it 


(e) cos n6 dB 


Combining equations (1) and (3) gives 


p+T! 

^i I si 
^n=-i - 

J -n 


cos n0 d0 


( 3 ) 


(4) 
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Since the integrand I is an even function of 6, where 1(6) = I(-0), 




2F. 

■7 


f 


(^: P) cos n0 d0 
9 - p 


( 5 ) 


These integrals are evaluated for both force components, F and F , at each node 
for all 25 values of n. Since 



( 6 ) 


the value of F^ in each case is the magnitude of the r or z force component assigned 
to that node. 


Centrifugal Body Force Representation 

There are two types of centrifugal body forces (inertia loads) on each element of 
volume in the ball. The first type of loading results from the rotation of the ball about 
its own spin axis at speed co. This rotational axis is the 2 - axis shown in figure 3, The 
second type of centrifugal load is produced by the orbital motion of the ball about the 
shaft centerline at speed For purposes of computation, the ball spin axis was con- 
sidered to be always parallel to the shaft axis. 

Spin loads . - Each ring element is centered about the ball z-axis and is subjected to 
radial centrifugal forces which depend on the material density p, cross-isectional area 
A, element radial position r from the z-axis, and the speed of rotation u) about the 
z-axis. For each ring element, the radially directed force per unit of circumference is 

"'spin ■ = 0 1 

which is directed as shown in figure 7. 

Orbital loads . - The orbital load per unit of circumference at angle 9, for a ring 
element located at ball radius r, and at distance R(0) from the shaft centerline, is of 
magnitude 


‘'orbital ' = P I 


( 8 ) 
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and is directed as shown in figure 7. Note that is dependent on e, since R(0) 

for any given r value can be described by 


R(0) = ycR + r cos 9)^ + (r sin 6)^ (9) 

where R is the bearing pitch radius (fig. 7). If r is small compared to R, R(e) is 
approximately equal to the pitch radius. 

The direction for also dependent on 9, but this load can be resolved 

Into two components, along the r and 9 directions, as shown in figure 8. 

The magnitudes of the two components can be approximated by 

® ''orbital * (10) 


'o * 'orbital « (11) 

This approximation assumes that the at^le (fig. 7) is small compared to 9. This 
approximation is valid when R is much greater than r. 

The importance of equations (10) and (11) is related to the Fourier representation of 
forces, which is used in this investigation. Note that if all coefficients except the first 
are zero in 


and 


then 


and 


F^(0) = 2^ \cosne+-^ 

n=l 


( 12 ) 


F,(0) 



cos n0 + 


T 


(13) 


Fj.(0) = bj cos 9 


(14) 


Fg(8) = dj sin 9 


(15) 
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Equations (14) and (15) match equations (10) and (11), respectively, if and d| are 
set equal to as given by equation (8). 

This method of approximating the orbital body loads is important since only Fourier 
coefficients are utilized by the computer prc^ram to represent loads. All orbital motion 
body loads are zero except for those in the set of elasticity equations for which n=l, 
whereas contact loads are nonzero for all values of n. 

In the actual calculations, the stresses due to the spin loads (since they do not de- 
pend on the angle 0) were calculated separately from those stresses due to the contact 
and orbital loads. The total stress then, by superposition, is equal to the sum of the two 
stress values for any given element. 


RESULTS AND DISCUSSION 

Stresses in the ball were computed for four cases, as stated previously, for a ball 
operating in a 125- millimeter -bore bearing. Two cases were for the regular drilled 
ball as tested in reference 8, and two cases were for a drilled ball with a 1. 5-millimeter- 
(0. 06-in. -) thick stiffening web. The actual force components applied to the balls in the 
analysis are given in table III. The material properties for ATSl M-50 steel utilized for 
the stress calculations are given in table IV. Note that with a ball radius of 10. 3 milli- 
meters (0. 406 in.) and a bearing pitch radius of 78.7 millimeters (3. 1 in.) the maximum 
angle (see fig. 7) is 7.5° (when 0 = 90°). Therefore, the approximate equations (10) 
and (11), based on the assumption that the angle is small compared to the angle 0, 
should be valid. 

The regular drilled-ball cases were computed first to determine if the values of cal- 
culated stresses were reasonable. Then the drilled-ball-with-web- cases were calculated 
to determine the effect on stress magnitude of the addition of a stiffening member at the 
center of the hole. The computer program output included the orthogonal stresses for 
each finite element at several values of the angle 9 (measured from the outer-race con- 
tact load position). These stresses are radial stress (r direction), transverse stress 
(z direction), and tangential stress (0 direction). All calculations were made in U. S. 
customary units and converted to SI units for the report. 


Drilled Ball 


The drilled-ball failures noted in reference 8 were attributed to flexural fatigue 
cracks that initiated in the bore of the ball. Therefore, the calculated stresses along the 
bore (at 0 = 0°) are plotted as shown in figure 9, for both the 11 000 and 22 000-newtons 
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(2500- and 5 000- lb) load conditions. 

The maximum bore stress is in the tangential ( 9 ) direction and occurs in element 115 
(see fig. 4) for the 11 000- newton (2500-lb) load condition. This element is located about 
1 millimeter (0. 04 in. ) from the centerline of the ball on the applied load side. For the 
22 000 newton (5000- lb) load condition, the maximum bore stress occurred in element 
147, which is about 3 millimeters (0. 1 in.) from the ball centerline. 

The values of tangential and transverse stress for these two elements with maximum 
bore stress are plotted as a function of the angle 6 in figure 10, Note the large revers- 
ing tangential stress that occurs for these elements as the ball rotates through one-half 
revolution. There are two stress reversals then for each complete ball revolution, one 
being larger than the other because of the difference in ball contact loads at the inner and 
outer race. 

Experimental work with hollow rollers in reference 9 resulted in a recommendation 
of about 500x10 newtons per square meter (70 000 psi) as the maximum permissible 
tangential stress level for bearings with rollers subjected to bending fatigue. Refer- 
ence 10 notes an endurance limit of 620x10^ newtons per square meter (90 000 psi) for 
a smooth AISI M-50 specimen. The maximum calculated tangential bore stress for the 
22 000-newton (5000-lb) case is about 470x10 nev/tons per square meter (70 000 psi) in 
tension, and the maximum value for the 11 000- newton (2500-lb) case is about 350x10 
newtons per square meter (50 000 psi). Note that, while the bearing load is doubled, the 
stress levels in the ball change by a factor of less than 1, 5. 

In bending fatigue, the chaj^e in stress (stress range) and the maximum effective 
stress are important in establishing an endurance limit for the material. For the 
22 000-newton (5000-lb) case, the largest total change in tangential bore stress is about 
620x10^ newtons per square meter (90 000 psi) and the maximum von Mises - Hencky 
effective stress at the bore is approximately 410x10 newtons per square meter (60 000 
psi). 

Since the bearings in reference 8 did operate for some time before failure, these 
calculated values for the stresses in the ball were judged to be reasonable. This judg- 
ment in turn added confidence to the method utilized for the stress calculation. Therefore, 
for comparison, the stresses for the ball with web were also calculated. 


Ball With Web 

The ball-with-web geometry was chosen as a design to reduce flexure in the drilled 
ball by adding a stiffening web at the center of the ball. Stress calculations were made to 
determine the effectiveness of this design change. 

The calculated tangential stress along the hole bore (and across the web) is plotted in 
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figure 11 for both the 11 000- and 22 000-newton (2500- and 5000- lb) conditions. The 
presence of the web not only reduced the stress level at the bore from that of the regular 
drilled ball, but also changed the character of the curves as well. The stress values 
shown are for 0= 0°. The maximum bore stress occurred in element 54 (fig. 5) for 
both load conditions. 

Since the changes in stress per revolution of the ball are important for flexural fa- 
tigue, values of tangential stress for element 54 were plotted as a function of the angle 
0 in figure 12. The stresses are shown for both the 11 000- and 22 000- newton (2500- 
and 5000- lb) conditions. 

The maximum calculated value of the tangential stress at the bore for the webbed 
ball is about 300x10^ newtons per square meter (45 000 psi) in tension for the 22 000- 
newton (5000-lb) case and about 240x10 newtons per square meter (35 000 psi) for the 
11 000-newton (2500-lb) case. The total change in stress with angular position at the 
22 000-newton (5000-lb) load is about 300x10 newtons per square meter (45 000 psi) for 
element 54. Note that this stress change for the webbed ball is only about one- half that 
calculated for the regular drilled ball (fig. 10). Also, note that the maximum tangential 
stress is only a little more than one-half that of the regular drilled ball. The maximum 
von Mises - Hencky effective stress at the bore (at 0 = 0°) was about 250x10® newtons 
per square meter (37 000 psi). Thus, the presence of the web significantly diminishes 
the large reversing stress that occurs for the drilled ball. 

For better visualization, a plot of the tangential stress at the bore, for both the 
drilled ball and the webbed ball, is presented on polar coordinates in figure 13 for the 
11 000-newton (2500-lb) condition. The stress values are for element 115 for the drilled 
ball (fig. 10) and element 54 for the webbed ball (fig. 12), Note that a tensile stress is 
plotted radially inward. Since strain is proportional to the stress, this plot gives a rep- 
resentation of the deformation of the ball at the bore surface due to the inner- and outer- 
race contact loads. 


SUMMARY OF RESULTS 

A finite- element technique was used to calculate the stresses in a drilled ball and in 
a ball with a stiffenir^ web. The ball with web was analyzed to determine the effect of a 
geometry change on the stress. The analysis was made for a 20. 6- millimeter- (0. 8125- 
in. -) diameter drilled ball, with a 12.6 millimeter- (0. 496-in.-) diameter hole (50-per- 
cent weight reduction). The web used was 1.5 millimeters (0.06 in.) thick, which re- 
sulted in a ball with a 46-percent weight reduction. The calculations were made for the 
conditions of a ball operating in a 125- millimeter-bore ball bearing, with a shaft speed 
of 24 000 rpm and bearing thrust loads of 11 000 and 22 000 newtons (2500 and 5000 lb). 
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The following results were obtained: 

1 . The maximum stress calculated for an element at the hole bore was in the tang- 
ential direction and occurred at the outer-race contact load position. At the 22 000- 
newton (5000 lb) bearing load condition, this maximum bore stress ranged from 470x10 
newtons per square meter (70 000 psi) in tension to 150x10® newtons per square meter 
(20 000 psi) in compression for one revolution of the drilled ball. 

2. The introduction of a stiffenii^ web significantly reduces the reversing tangential 
stress at the hole bore. The total stress change per revolution for the webbed ball was 
about one-half that for the regular drilled ball. 

3. The analysis, made for a bearir^ ball with nonaxisymmetric loading, resulted in 
reasonable bore stress values. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 15, 1973, 

501-24. 
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TABLE L ' BEARING SPECIFICATIONS 


Bearing outside diameter, mm (in,) 190 (7.48) 

Bearing inside diameter, mm (in.) 125 (4.92) 

Bearing width, mm (in.) 52.5 (1.28) 

Outer -race curvature 0.52i 

Inner-race curvature 0,52 

Number of balls 21 

Ball diameter, mm (in. ) 20.6 (0.8125) 

Ball hole diameter, mm (in.) ■ 12.6 (0.496) 

Bearing pitch diameter, mm (in.) 157.5(6,2) 

Retainer type piece, with stub pins 

Retainer material AISI 4340, silver plated 

Race and ball material AISI CVM-50 






































TABLE IV. - MATERIAL PROPERTIES OF AISI M-50 STEEL USED IN CALCULATIONS 


Density, kg/m^ (Ib/in.^) 8.0x10^(0.29) 

Young's modulus, N/m^ (psi) 2.0x10^^ (29x10®) 

Poisson’s ratio 0. 25 



Rgure 1. - Drilled-ball bearing with 125-millimeter-diameter bore tested in refer- 
ence 8 and used for analysis in present report. 



(a) Drilled ball. 

Figure 2. - Cross section of ball 


(b) Ball with web. 

configurations used for stress analysis. 
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Figure 3, - Ball contact loads and angles at operating conditions and cylindrical coordinate 
system used for analysis. 



figure 4. - Finite-element mesh for drilled ball in r, Z-plane. 









Elemental tensile 

stress, psi Elemental tensile stress, psi 


-Tangential stress 


Bearing toad. Element 
N m 

11 000 ( 2500) 115 {fig. 4) 

22 000 (5000) 147 (fig. 4) 


V 

300 — \ 


200H \\ 


V Transverse stress 


0 20 40 60 80 100 120 14D m 18i 

Angle from outer -race load position, 0, deg 

Figure 10. - Stress at bore of drilled ball as function of angle from outer-race load position 9. Elements for maxi 
mum stress at 0 =* 0^. 
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Figure 11. - Tangential stress along bore of ball with web at angle from outer-race load position 0 of 0° for two 





Figure 12. - Tangential stress at bore ot ball with web, as function of angle from outer- race load position 0 for 
element with maximum bore stress at 0 = oP {element 54, fig. 5). 
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Figure 13. - Tangential stress at bore of bearing ball as function of angle from outer- race load 
position 6. Bearing load, 11 OOO newtons (2500 lb). 
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